Entropic droplets and activated events near the glass transition of 

a random heteropolymer 

Shoji Takada and Peter G. Wolynes 
School of Chemical Sciences, University of Illinois, Urbana, IL, 61801 

(February 1, 2008) 

Abstract 

The barriers between metastable states near the glass transition of a random 
heteropolymer are studied using replicas by describing inhomogeneous states. 
The instanton solution for a replica space free energy functional is found 
numerically to estimate the size of activation barriers and of the critical nuclei 
themselves between the dynamic and the static glass transition temperatures. 
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Quantifying escape from configurational traps on a rugged free energy landscape is im- 
portant for understanding the dynamics of spin glasses structural glasses 0, and folded 
proteins 0, and for protein folding |Q. A first step towards understanding barriers is to 
appreciate the organization of the stable minima This has only been carried out com- 
pletely for infinite range spin glasses. That organizational structure inspires many of the 
dynamical theories Barriers in mean field spin glasses scale with the system size, 

but the finite range of interactions allows escape from traps through localized reconfigura- 
tions with finite barriers in the thermodynamic limit. The Vogel-Fulcher law for viscosity 



of structural glasses has been explained through such a mechanism |jTO|. Various mean field 
theories of structural glasses resemble those for spin glasses lacking reflection symmetry. In 
strict mean field theory these models undergo a dynamical transition at a temperature Ta, 
where a macroscopic number of frozen free energy minima appear, and a static transition 
at Tk, the Kauzmann temperature, where the configurational entropy of the minima disap- 
pears. Kirkpatrick and Wolynes (KW) pointed out that that individual free energy minima 
between Ta and Tk will be inherently unstable for short range interaction models because 
of entropic droplets. The extensive configurational entropy provides a driving force for a 
localized region in a local minimum to reconfigure and escape the trap. Their analysis gave 
a modified Vogel-Fulcher law while a later scaling picture incorporating entropic droplets 
gave precisely the usual form used empirically [0. Parisi has presented a novel instanton 



argument in replica space yielding the original KW form . Here, we use replica instanton 



calculations to quantify reconfiguration barriers for the random heteropolymer. 

Reconfiguration barriers determine the configurational diffusion coefficients that enter 
the theory of protein folding times [1,0]. At low T, trap escape is rate limiting. Whether 



the escape barrier is extensive can be tested experimentally in a crude fashion since pro- 
teins are mesoscopic with finite chain length A^. The size scaling is presently controversial 
TSHTSl. The barriers computed in the strict mean field limit of uniform transitions are 



of order ~ O.lA^fc^Tx at the static transition Tk using parameters fit to the lattice model 
thermodynamics This seems to be consistent with recent simulation results at low tem- 



perature [|T4|. Our calculations suggest entropic droplets are often large in the random 
heteropolymer so that the mean field arguments are a good starting point for mesoscopic 
systems of the size of the smaller proteins. One caveat is that polymers have additional 
entanglement constraints not present in structural or spin glasses. These are neglected here. 

We focus on the temperature range near the glass-like transition of a finite-size het- 
eropolymer in a poor solvent so that it is collapsed. The droplet we deal with is a small 
globular region of the polymer that may take on many configurations, local in space but not 
necessarily in sequence, buried in a remaining frozen glassy portion. We utilize replica for- 
malism and derive a Landau-like free energy functional in terms of the Debye- Waller factor 
for a residue which plays the role of an Edwards- Anderson order parameter which is taken 
as spatially varying. 

We introduce the standard bead contact Hamiltonian for a random heteropolymer which 
includes finite range random interactions between monomers; 

H = ksTY. ^''^\~/^' + iT^b^^e-^'^-'^^'/^' + I E e-[(--^^)^-^(r--)1/-^ (1) 

where represents positions of Kuhn segments = 1 ~ A^), a is the Kuhn length, 
6jj and c are the second and third virial coefficients, respectively, and a is the char- 
acteristic length of interactions. 6jj are Gaussian random variables with distribution, 
P(6,,) = (27r62)-i/2exp \-{h,^ -h,f I2h\ 

The free energy is first averaged over this random interaction using stan- 
dard replica formalism [^], to obtain an effective Hamiltonian, ifefr = -f^o + -f^i + 
H2 where = fc^T X^ia (r^^+i — r°)^/2a^ represents the elasticity. Hi = l/2(6o — 
/?6V2)Ei^,aexp[-(r--r-)Va2] + c/6E^^i^fcae-[('■"-'•?)'+^"-'■^)']/'^' are effective ho- 
mopolymeric interactions and n is the number of replicas. The inter-replica interaction 
part H2, responsible for breaking ergodicity, can be written in terms of an overlap order 
parameter function Qap defined by Qapi^'i, ^2) = J2i (^(fi — r?)<5(r2 — rf ) 0]) where a and 
are replica indices; 



3 



dridr2dr3dr4 



-[(ri-r3)2+(r2-r4)2]/a2 



(2) 



Exact integration over the bead variables is difficult, but a variational approach extended 
into replica space |T^|8|,P| can be used with a reference Hamiltonian ifj-ef to calculate a 



variational free energy Fvar by finding extrema of Fvar = —ksTlnZj-ef + (i^eff — -f^ref)- ^ref 
is the partition function for ifj-cf and (■ ■ ■) means the average with ifrcf- The physical free 
energy is \imn-,o F*^^/n, where * means its extreme value. We use the same reference 
Hamiltonian as ref.I (with C = in ref.I), 

f^H^ef = E(rr+i - rtr/{2a') + B Y^ir?? + D d^p{vt - rf )^ (3) 

where B measures the confinement to a globule and D, and da/s are variational parameters 
specifying the vibrational freedom in a minimum and the replica symmetry breaking related 
to the configurational entropy, respectively. We assume dap has the same structure as mean 
field Potts spin glasses; n replicas are divided into n/m groups, each of which has size 
m and the matrix element d^/s is 1 if a and f3 {a ^ (3) belong to the same group and 
otherwise. It is straightforward, though cumbersome, to obtain Fvar as a function of these 
parameters. As in Ref.I the homogeneous glassy state characterized by a large constant D 
(i.e. D ^ A ^ B) in the reference Hamiltonian yields an asymptotic high D expression 
for the free energy of the homogeneous glassy phase. A globally appropriate homogeneous 
expression can be obtained by interpolating to the D = limit yielding the non-gradient 
terms in eq.(5) (see below). 

An entropic droplet is described by an inhomogeneous situation where part of the polymer 
is trapped in a particular metastable state while another part can be in any minimum. The 
exterior of the droplet has large D describing the part trapped in a particular metastable 
region while the interior with small D represents a region which can exist in multiple states. 
For explicitness, we shall assume the interaction range a is small compared to the scale on 
which D varies. Expanding the quadratic order parameter interaction in eq.(^ around the 



midpoint of the two coordinates yields a (V-D)^ term describing the surface tension between 
mean field minima |T^. Defining a dimensionless Debye- Waller factor y = 2mDcr^/2, the 



free energy functional becomes in this approximation F[y(r)] = -FGiobuie + / dr/ {y (r)), where 
-^Globule is the constant free energy for the globule and 
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/ (y (r)) = TipksT In 
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(m - 1)— 



^y + lj 128 {y + 1)7/2 

(4) 



is the inhomogeneous free energy density. p(r) is the density of beads and u = (7r(T^/2)^/^. 
The surface tension is proportional to b"^ and originates from the randomness of the inter- 
actions. Although not strictly true for the Gaussian confinment model, we have ignored 
derivatives of p, which are expected to small in a well collapsed state with excluded volume. 
In the functional, m y and p are functions of r. The former two play the main roles in glassy 
behavior. 

The uniform solutions are essentially the same as in ref.I for a contact interaction model. 
Since p does not exhibit any peculiar behavior, we use the so-called volume approximation 
where the monomer density pQ,(r) = p is taken inside polymer and zero otherwise. Putting 
y{r) = y, we first maximize with respect to m for each y (denoting it as m*{y)). Fig.|I| shows 
f{Tn*{y),y) for several temperatures. In the high temperature limit, there is no saddle 
solution and only the globule state y = is stable. The dynamic glass transition takes place 
at the temperature Ta where a non-zero saddle solution y > appears. This describes a 
glassy trapped state. At temperatures below Ta, there are two minima in f{y), one for the 
globule {y = 0), which is thermodynamically equivalent to the sum over metastable minima, 
and the other for a particular metastable glassy minimum, {y = yQ > 0). As temperature 
decreases, the free energies for the two solutions approach each other and become degenerate 



at the Kauzmann temperature Tk where the static transition occurs. As in spin glasses []T8 
the uniform solution yields an estimate of a lower bound for the barriers between two lowest 
minima as the local maximum of f {m* (y) , y) , which we denote as occurring at y^. Starting 
from zero at Ta the uniform solution barrier height is extensive and grows as temperature 



decreases and saturates around Tk 0. T^/T^ evaluated with typical parameters for flexible 
polymers {{2a/ a)^ = 4, /w = 1) is about 1.4. 

The globule state between Ta and Tk obtained by the ordinary replica theory represents 
a weighted sum of multiple local minima. Thus, the difference between the most probable 
free energy of a minimum and the globule free energy corresponds to the logarithm of the 
degeneracy of local minima, i.e., an extensive basin configurational entropy. For the short 
range infinite system between Ta and Tk local minima are not separated by a thermodynam- 
ically large barrier but a finite one. KW argue the free energy of an entropic droplet with a 
radius i? is T ~ Si?^ — TscR^, where S is the surface tension and Sc is the configurational 
entropy density. Since Sc disappears at Tk, the critical size of the droplet diverges at Tk 
where the static glass transition occurs. 

Inhomogeneous saddle point solutions of the free energy functional in the replica formu- 
lation give a microscopic theory of the barrier analogous to the KW result for Ta > T > Tk. 
The random interactions are integrated out and so the interface of the replicated droplet 
will be spherical, minimizing the surface energy. In spherical coordinates, the stationary 
phase condition for y{r) leads to an Euler-Lagrange equation [^, that corresponds to the 
Newton equation of a dissipative system, where time is r and a coordinate is y. For an 
infinite polymer, the appropriate solution has the boundary conditions dy/dr\r=o = and 
y{oo) = yQ. Since the dissipation diverges at r = 0, the solution stays at y = at the begin- 
ning of trajectory (r ~ 0), falls after a while with an infinitesimal initial velocity, and stops 
aX y = yQ. Since the dissipation monotonically decreases with r, there is one such trajectory 
for each temperature. The critical radius r^ defined by the value of r where ?/(rc) = is 
shown as a function of temperature in Fig.^ shows the critical free energy of the 

droplet in an infinite system as well as that for a uniform transition for a finite but large 
A^. In these figures, we indicate by dashed lines the results obtained by the so-called thin 



wall approximation [|T9| , which assumes a plane domain wall giving a free energy of the KW 
form. 

For protein-like parameters, just below Ta, the free energy barrier for a droplet in an infi- 
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nite system is higher than the barrier for a uniform transition of a polymer with a moderate 
size (e.g. 1000- mer). The interface for the droplet of small radius costs a considerable surface 
energy reflecting the difficulty of healing the interface. A polymer of this size escapes from 
traps through a nearly uniform reconfiguration. Decreasing temperature, a localized escape 
mechanism begins to dominate (Fig.^ for a sufficiently large polymer, but the critical ra- 
dius then increases upon further cooling until it again reaches the whole size of the polymer, 
because the entropic driving force for transition finally disappears at Tk. Thus, very close 
to Tk, the spatially uniform transition again dominates. The width of the spherical domain 
wall is ~ 3cr. This is not much smaller than the radius of the smaller lattice models for 
proteins. The thin wall approximation is accurate for larger droplets (lower temperature), 
but it drastically underestimates the free energy for small droplets. The numerical accuracy 
of these results, of course, depends on how well the real system is approximated by the 
contact model but the trends should be robust. 

We have described droplets that are local in space. There is an alternative escape route; 
an entropic droplet local in sequence. This can be treated using a similar reference Hamil- 
tonian to eq.(|^), but with Di depending on sequence number i. Straightforward calculation 
shows that the activation barrier for the sequentially localized droplet is proportional to the 
size of system N between Ta and Tk, while for droplets local in space the activation energy 
is independent of except very close to Tk. Entropic droplets local in sequence therefore 
do not change the story much. Droplets local both in sequence and space correspond for 
glassy traps with the foldons of a minimally frustrated system and may be relevant when 
topological constraints are considered ||20[ . 



The explicit droplet solutions discussed here are oversimplified. For structural glass and 
Potts glasses the interface has a more complex structure ||TT|. Wetting due to multiple 
states in the interface reduces the effective surface tension significantly. This leads to the 
Vogel-Fulcher dependence AT ~ (T — Tk)^^ We have not yet succeeded in quantifying 
this wetting phenomenon in the replica instanton formalism. Inhomogeneous states of a 
polymer (and short range Potts- type models, in general) may not, rigorously speaking, be 



described by the simple one level RSB scheme but require a P{q) with finite width peaks. A 
treatment of the interface like that for the short range SK model [|12| may incorporate the 
wetting effect. 

We are grateful to Jin Wang and John Portman for useful discussions. S.T. is a Postdoc- 
toral Fellow for Research Abroad of the Japan Society for the Promotion of Science. P.G.W. 
is supported by NIH grant PHS 1 ROl GM44557. 
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FIGURES 

FIG. 1. Free energy density for the uniform solution as a function of y = ImDa^ jl for several 
temperatures. The value of m is maximized for each y. Parameters used are {2a /aY' = 4 and 
pv = 1 and Temperatures are T = 0.326, 0.2926 ~ Ta, 0.246 and 0.20396 ~ Tk. 

FIG. 2. Critical radii of droplet plotted with respect to T. The solid (dashed) curve represents 
the spherical droplet (that by the thin approximation). Parameters used are the same as Fig.||. 

FIG. 3. Free energy barriers for droplets and uniform transitions with respect to T. Each curve 
is for l)the droplet, 2)the droplet by the thin approximation, 3)the 100-mer uniform transition and 
4)that of 1000-mer. Parameters used are the same as Fig.||. 
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